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Abstract 

It has recently been suggested that collisions of transversely polarized protons and 
antiprotons at the GSI could be used to determine the nucleon’s transversity densities 
from measurements of the double-spin asymmetry for the Drell-Yan process. We 
analyze the role of higher-order perturbative QCD corrections in this kinematic regime, 
in terms of the available hxed-order contributions as well as of all-order soft-gluon 
resummations. We hnd that the combined perturbative corrections to the individual 
unpolarized and transversely polarized cross sections are large. We trace these large 
enhancements to soft gluon emission near partonic threshold, and we suggest that 
with a physically-motivated cut-off enhancements beyond lowest order are moderated 
relative to resummed perturbation theory, but still signihcant. The unpolarized 
dilepton cross section for the GSI kinematics may therefore provide information on 
the relation of perturbative and nonperturbative dynamics in hadronic scattering. The 
spin asymmetry turns out to be rather robust, relatively insensitive to higher orders, 
resummation, and the cut-offs. 


1 Introduction 


The partonic structure of polarized nucleons at the leading-twist level is characterized by the 
unpolarized, longitudinally polarized, and transversely polarized parton distribution functions /, 
A/, and 6f, respectively p. In contrast to the distributions / and A/, we have essentially no 
knowledge from experiment so far about the transversity distributions 6f, even though there are 
now first indications j2] that some of them are non-vanishing. The 6f were first introduced in p. 
They are dehned as PPUP the difference of probabilities for finding a parton of flavor / at 
scale /i and light-cone momentum fraction x with its spin aligned (|t) or anti-aligned (jt) to that 
of the transversely polarized nucleon: 

(5/(a:,/i) =/lt(a;,/i) -. (1) 


By virtue of factorization theorems inip, the parton densities can be probed universally in a 
variety of inelastic scattering processes for which it is possible to separate (“factorize”) the long¬ 
distance physics relating to nucleon structure from a partonic short-distance scattering that can 
be calculated in QCD perturbation theory. It was realized a long time ago pinna that due to its 
chirally-odd structure, transversity decouples from inclusive deeply-inelastic scattering, but that 
inelastic collisions of two transversely polarized nucleons should offer good possibilities to access 
transversity. In particular, the Drell-Yan processes pp —>■ pp —>■ l^l~X {I = e,fi) were 

identified as promising sources of information on transversity iHiiniiini. This is so because there 
is no gluon transversity distribution at leading twist pp. For the Drell-Yan process, the lowest- 
order partonic process is qq ^ J*, with gluonic contributions to the unpolarized cross section in 
the denominator of the transverse double-spin asymmetry 


_ gTT _ 

(jTT -|- Q-U 


( 2 ) 


only arising as higher-order corrections. Therefore, A^t may be sizable for the Drell-Yan process, 
in contrast to other hadronic processes such as high-transverse-momentum prompt photon and 
jet production p cn ini cni [H] which in the unpolarized case are largely driven by gluons in the 
initial state and are hence expected to have a very suppressed Att- 


Clean information on transversity should be gathered from polarized proton-proton collisions at 
the BNL Relativistic Heavy Ion Collider (RHIC) where the Drell-Yan process is a major focus jl5j . 
In pp collisions, however, the Drell-Yan process probes products of valence quark and sea antiquark 
distributions. It is possible that antiquarks in the nucleon carry only little transverse polarization 
since, due to the absence of a gluon transversity, a source for the perturbative generation of 
transversity sea quarks from g ^ qq splitting is missing. In addition, at RHIC energies and for 
Drell-Yan masses of a few GeV, the partonic momentum fractions are fairly small, so that the 
denominator of A^t is large due to the small-a; rise of the unpolarized sea quark distributions. 
Thus, even for the Drell-Yan process, the spin asymmetry A-tt at RHIC will probably be at most 
a few per cent, as theoretical studies have shown pcni unl¬ 
it has recently been proposed to add polarization to planned pp collision experiments at the 
GSI, and to perform measurements of Att for the Drell-Yan process dUlIHlIinilini- This is a 
very exciting idea, since unique information on transversity could be obtained in this way. The 
results would be complementary to what can be obtained from RHIG measurements. First of all. 
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in pp collisions the Drell-Yan process mainly probes products of two quark densities, 6q x 6q, since 
the distribution of antiquarks in antiprotons equals that of quarks in the proton. In addition, 
kinematics in the proposed experiments are such that rather large partonic momentum fractions, 
X ~ 0.5, are probed. One therefore accesses the valence region of the nucleon. Estimates dniEoiEi] 
for the GSI PAX and ASSIA experiments show that the expected spin asymmetry Att should be 
very large, of order 40% or more. 


It is important, however, to keep in mind that the kinematic region to be accessed in the hrst 
stage of the GSI experiments, with Drell-Yan masses M of 1 — 4 GeV or so, but a center-of-mass 
energy of only a/ 5 ~ 5.3 GeV for the baseline hxed-target program, is not really the “classic” 
regime where parton model ideas, factorization, and perturbative QGD are a priori expected to 
provide adequate descriptions. This is of course crucial since the interpretation of A^t in terms 
of transversity relies exactly on these concepts. To be more precise, at high energies and large 
dilepton invariant mass M the cross section factorizes inmi into convolutions of parton densities 
and perturbative partonic hard-scattering cross sections, as mentioned above. Schematically, 




da 


fa® fb<^ 

a^b 


M'^ddab 

dM^ 



p 


( 3 ) 


For simplicity, we have considered here the unpolarized cross section, and we have also integrated 
over the rapidity of the lepton pair and only focused on the total Drell-Yan cross section. We 
also have not written out the precise form of the convolutions, which will be given below. For the 
moment, we are only interested in the important features visible in Eq. dBD- The quantities one 
wants to determine from measurement of the left-hand-side of Eq. (H are the parton distributions 
/a, fb- The partonic cross sections, dab, for the reactions ab —>• 7*X may be calculated in QGD 
perturbation theory. Their expansion in terms of the strong coupling constant as{M) reads 


ddab = dd^°J + 


as{M) 


TT 


A” + 




dd 


( 2 ) 

ab 


+ . . . , 


( 4 ) 


corresponding to lowest order (LO), next-to-leading order (NLO), and so forth. The earlier studies 
113120] for the Drell-Yan process at GSI energies used LO hard-scattering cross sections to estimate 
the expected spin asymmetries. Depending on kinematics, however, the higher-order corrections 
may be very important. As we will show below, this is the case for the planned GSI measurements. 
In addition, as indicated in Eq. m, factorization of the hadronic cross section in terms of twist-2 
distributions is of course not exact, but holds only to leading power in M. There are corrections 
to the (dimensionless) cross section M'^da/dM‘^ that are down by inverse powers of the hard scale, 
that is, of the form {X/MY with some p and some hadronic mass scale A [221 • These power 
corrections will also depend on r = M'^/S and are generally expected to increase with increasing 
r. The measured spin asymmetry A^t can only be reliably interpreted in terms of the transversity 
densities if the higher order and power corrections can either be shown to be small in the accessible 
kinematic domain, and/or if they are sufficiently well understood. The aim of this paper is to 
address primarily the question of how large the higher-order QGD corrections to the Drell-Yan 
process are in the GSI kinematic regime. To this end, we will apply the technique of threshold 
resummation, to which we now turn. 

As we mentioned above, r is typically very large for the GSI kinematics, 0.2 / r < 0.7. This is 
a region where higher-order corrections to the partonic cross sections are particularly important. 
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r = 1 sets a threshold for the reaction, and as r increases toward unity, very little phase space 
for real gluon radiation remains in the partonic process, since most of the initial partonic energy 
is used to produce the virtual photon. Virtual and real-emission diagrams then become strongly 
imbalanced, and the infrared cancellations leave behind large logarithmic higher order corrections 
to the partonic cross sections, the so-called threshold logarithms. At the kth order in perturbation 
theory, the leading logarithms are of the form ln^^“^(l — z)/{l — z), where z = r/xaXh is the 
partonic analogue of r. For sufficiently large z, perturbative calculations to hxed order in 
become unreliable, since the double logarithms compensate the smallness of as{M) even if M is 
of the order of a few GeV. The fact that the parton distributions are steeply falling functions 
of the momentum fractions Xafi means that the threshold region is actually emphasized in the 
cross section, even if r itself is still rather far away from one. If r is close to unity, as is the 
case for much of the GSl kinematics, the region of large z < 1 completely dominates, and it is 
crucial that the terms ln^^“^(l — z)/{l — z) be resummed to all orders in ag. Such a “threshold 
resummation” was originally developed for the Drell-Yan process and subsequently applied 

to a variety of more involved partonic processes in QGD h It turns out that the soft-gluon effects 
exponentiate, not in z-space directly, but in Mellin-V moment space, where N is the Mellin 
moment conjugate to z. The leading logarithms (LL) in the exponent are of the form ln^'''^(V), 
subleading logarithms (next-to-leading logarithms (NLL)) of Q;^ln^(V). In this paper we will use 
NLL resummed perturbation theory to analyze the importance of higher-order corrections to the 
Drell-Yan cross section in the kinematic regime to be explored at the GSl. 

There is a close relation between resummation and the nonperturbative power corrections. 
Taking into account nonleading logarithms and the running of the coupling, resummation always 
leads to a perturbative expression in which the scale of the coupling reflects the value of the 
transform variable. Because of the singularity of the perturbative effective coupling at Aqcd, the 
resulting expressions are ill-dehned j2H|- The analysis of these ambiguities for Drell-Yan cross 
sections j2niEn] suggests a series of nonperturbative corrections ra, generically suppressed by 
even powers of the pair mass M, but enhanced by the moment variable N, As we 

shall see, for GSl hxed-target energies, the effective values of N are so large that the hrst few 
power corrections will not suffice for the Drell-Yan cross section. We will therefore rely on a 
somewhat different approach, to be presented in more detail elsewhere, and cut off unphysical 
dependence on low momentum scales. The result for a large portion of dilepton masses M will be 
a cross section with moderated, but still signihcant enhancements relative even to next-to-next- 
to-leading order calculations, which we take as a conservative prediction based on perturbation 
theory. Experimental results on these cross sections should shed light on the interrelations between 
hxed-order, all-order and nonperturbative corrections in hadronic scattering. 

The remainder of this paper is organized as follows. Section |21 will present the basic framework 
for our calculations and will introduce the partonic threshold region. In sec. El we provide all 
ingredients for the NLL resummation of the threshold logarithms, and we also propose a new 
infrared-regulated expression for the form of nonperturbative corrections suggested by perturbative 
resummation. In sec. 0] we then present phenomenological results for the regions of interest in GSl 
measurements. 

^See, for example, the review in Ref. HH], and interesting recent work [2^1 that rederives some of these results 
in the context of the soft-collinear effective theory 
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2 Perturbative cross section and the threshold region 


The spin-dependent cross section for dilepton production by two transversely polarized hadrons 
is defined as 

rdda _ 1 / rda^^ r(icr^^\ , . 

drdcj) 2 \ drdcf) drdcf) ) ’ 

where the superscript || (||) denotes parallel (antiparallel) setting of the transverse spins of the 
incoming hadrons. We have used the customary Drell-Yan scaling variable r = M‘^/S with M 
the invariant mass of the lepton pair and S the center-of-mass energy squared. 0 is the azimuthal 
angle of one of the leptons, counted relative to the axis defined by the transverse polarizations. 
For simplicity, we have integrated over all rapidities of the lepton pair. At high M, the cross 
section factorizes into convolutions of the transversity distributions 6f with the corresponding 
transversely polarized partonic hard-scattering cross sections [ 7 ]: 


= f 

n J T 

a,b 


dx 


dxh 


- 6fa{Xa, 11 ) dfb{Xb, II ) duJah Z = xjXaXb, 

Jrixa V 


M 2 


TdSa{T) 
drdcj) 

u.,u 

( 6 ) 

where /i collectively denotes the factorization and renormalization scales, and where we will specify 
the (0) below. Due to the odd chirality of transversity, and since there is no transversity gluon 
density, qq annihilation is the only partonic channel, up to next-to-leading order (NLO)h Its cross 
section is calculated in QCD perturbation theory as a series in 


6uJab {z, r, as) = Su'fJ (z) + (z, r) + ) duj^f (z, r) 

TT V TT / 


+ 


(7) 


where r = M^/Since we will not consider very high energies and dilepton masses, only photons 
contribute as intermediate particles. The lowest-order (LO, 0{a^)) process thus is q^q^ —> 7 * —>■ 
for which 


« 2 e 2 


(^) = ^(1 - ^) , ( 0 ) = cos( 20 ) . 


( 8 ) 


The hrst-order term <50;^^ is known and reads Bn 

Tn(l — 2 ;)\ 2z\nz 3z\n^ z 


S^^qq{z,r) = Cb 


4z 


1-z 1- 

, -,' 5 ( 1 -^) +(tt^ 

3 y 


+ 2(l-z) 


^-0 


-I- -^(1 — z) ) Inr 


(9) 


where Ci? = 4/3 and the “-I-”-distribution is defined as 


[ dz [g{z)]^ f{z) = j dzg{z) {f{z) - /(!)) 
'0 Jo 


( 10 ) 


Since resummation is performed in Mellin moment space, we take a Mellin transform of the 
hadronic cross section: 

d5a^ _ ^ N-i^d5a 

dcj) Jq drdcj) 


( 11 ) 


^Starting from next-to-next-to-leading order (NNLO) there are contributions from qq scattering as well. 
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The cross section algebraically factorizes under moments, 


d6a 


N 


d(j) 








TT 


TT 




( 12 ) 


where the Mellin moments of the transversity distributions dq and the higher-order corrections 
duj^qq are dehned as usual, 


= / dxx^'^ ^5fq{x,ii^), 


N-l- 


dujf^’^ir) = / dzz^ ^5uj^q^{z,r) . 


(13) 


The moments of the hrst-order correction duj^qq in Eq. © read in the MS scheme EB: 


fcg’’")’-) = C, 


_iV(iV + 1) 

The sums appearing here are dehned by 


+ 2Sl{N) + 6 {S^{N) - C(3)) - 4 + -TT^ + - - 2S^{N) Inr 


N 


sm - E-i 


(14) 


(15) 


i=i 


Their analytic continuations to arbitrary Mellin-are 


S,(Af) = V'(A' + 1)+7 e, 

SsiN) = li/>"(iV + 1) + C(3) , (16) 

where 'ijj^z) is the digamma function, = 0.5772... is the Euler constant, and C(3) ~ 1.202057. 


We mention that formulas analogous to the above hold for the unpolarized case. The main 
difference is that in the unpolarized case beyond LO there are contributions from initial-state 
gluons to the Drell-Yan cross section. In the kinematic region we are interested in here, these are 
rather unimportant. All details for the unpolarized case to NLO may for example be found in 

Ref. [g. 


The threshold region corresponds to —*■ 1 or Y —>■ oo. At large N, the moments of the NLO 
correction become 


fo;y'"(r) = cp 


21„2(iv)-4 + L=+ (5 


2 ln(7V)l In r 



(17) 


where 

N = . (18) 

One can see the double-logarithmic corrections cx Q!sln^(iV) near threshold, associated with the 
logarithmic term r\j ln(l — 2 ;)/(l — 2 ;) in Eq. Q. At higher orders, there are corrections of the form 
af'ln^(W) with I < 2k. We emphasize that the behavior of the unpolarized partonic cross section 
near threshold is exactly the same as Eq. (HZD. This is related to the fact that the large logarithms 
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are due to the emission of soft gluons, which is spin-independent. Thanks to the simple structure 
of the LO Drell-Yan process, the constant (Y-independent) pieces which are partly associated 
with virtual corrections are identical as well for the transversely polarized and unpolarized cases, 
provided both are treated in the same factorization scheme. 

We now turn to the resummation of the leading and next-to-leading threshold logarithms to 
all orders in 


3 Resummed cross section 


In Mellin-moment space, threshold resummation for the Drell-Yan process results in the exponen¬ 
tiation of the soft-gluon corrections. To NLL the resummed formula is given in the MS scheme 
by 


Suqf {r, asifi)) = exp [Cq{r, asi/J,))] exp I 2 j dz ^ ^ J ^ -^Aq{asikT))> , (19) 


where 

+ 

TT ^ 

(v)^Y 

+... 


(20) 

with jUj: 

= Cf . = 1 Cf 

c. (8^ - 

VIS 

y) 

^Nf , 

9 ^ ’ 

(21) 


where Nf is the number of flavors and Ca = 3. The coefficient Cg(r, Q!s(/^)) collects mostly hard 
virtual corrections. It is a perturbative series and reads 


Cq{r,as{^J)) = (^-4 -F ^ ^ Inr^ + 0{al) . (22) 


We note that it was shown in j21] that these corrections also exponentiate. 

Eq. (jl9|l as it stands is ill-dehned because of the divergence in the perturbative running coupling 
as(fcr) at kr = Aqcd- The perturbative expansion of the expression shows factorial divergence, 
which in QCD corresponds to a power-like ambiguity of the series. It turns out, however, that 
the factorial divergence appears only at nonleading powers of momentum transfer. The large 
logarithms we are resumming arise in the region [ 21 ] z <1 — 1/N m the integrand in Eq. ()19|1 . 
One therefore hnds that to NLL they are contained in the simpler expression 






'M2/7V2 


A / /I \\A NkT 

A,{a.{kr))\n — 


+ 2 


dk^T 


'M 2 


^Ag(a^(fcr))lnY 

iVrp 


(23) 


for the second exponent in (uni). This form, to which we will return below, is used for “minimal” 
expansions jUj of the resummed exponent. 

^We note that the threshold resummation for the Drell-Yan process has been worked out even to next-to-next- 
to-leading logarithmic accuracy, see m 
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3.1 Exponents at NLL 


In the exponents, the large logarithms in N now occur only as single logarithms, of the form 
ln*^^^(iV) for the leading terms. Subleading terms are down by one or more powers of In(iV). 
Knowledge of the coefficients in Eq. (HI is enough to resum the full towers of LL terms 

ln^'''^(iV), and NLL ones aJln^(iV) in the exponent. With the coefficient Cq one then gains 
control of three towers of logarithms in the cross section, ln^*^(iV), ln^^~^(iV), aj ln^*^“^(iV). 

We now give the explicit formula for the expansion of the resummed exponent to NLL accuracy. 
From Eqs. (Hi,(El one hnds [21 EH] 

\i\5u/qA^{r,as{^i)) = Cq(r, as(/i)) + 21niV + 2h(^)(A,r) , (24) 


where 


The functions are given by 


A = bottsi^) IniV . 








^ [2A + (1 - 2A) ln(l - 2A)] , 


2'K‘^hQ 


[2A + ln(l - 2A)] + 


27r6|] |_ 


2A + ln(l - 2A) + - ln2(l - 2A) 


4 ( 1 ) (lA 

—— [2A + ln(l — 2A)] ln(r)- - — - - In N ln(r) , 

27ron TT 


where 


6 n = 


hi = 


127r 

1 


247r2 


{lieA - 2Nf) , 
{l7Cl-5CANf-SCpNf) . 


(25) 

(26) 


(27) 


(28) 


The function contains all LL terms in the perturbative series, while is of NLL only. 
We note that the resummed exponent depends on the factorization scales in such a way that 
it will compensate the scale dependence (evolution) of the parton distributions. This feature is 
represented by the last term in (HI). One therefore expects a decrease in scale dependence of the 
cross section from resummation. The remaining p-dependence in the second to last term in El 
results from writing the strong coupling constant as 


asikr) 


(h) 

1 + boasifi) ln(/c|,//i2) 


+ O {asififiasin) ln(/c|/p^))”) 


(29) 


when doing the NLL expansion of the exponent. For this term, /i represents the renormalization 
scale. 


As was shown in Refs. El EH, it is possible to improve the above formula slightly and to also 
correctly take into account certain subleading terms in the resummation. To this end, we rewrite 
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Eqs. (HM-TITl as 


\n5ulf^{r,as{fJ.)) 


— [2A + ln(l - 2A)] 

71 bo 




+ 


71 




+ [-2AWlniV-5W] f 


TT 


Tlbn 


where = —“iCpI^. The last term in Eq. (IHOD is the LL expansion of the term 


Since the factor 


-2^^“’ IniV - bJ‘> 


(30) 


(31) 


is the large-iV limit of the moments of the LO DGLAP 

splitting function for transversity, the term in Eq. m may be viewed as an evolution of 

the parton distributions between scales /i and M/N. This suggests to modify the resummation 
by replacing IHOIIHHI 


[- 2 AlW InAT-^W] 


6P^ = C 


-25i(iV) 


(32) 


in Eq. PI|). To see the improvement resulting from this, we expand the resummed formula in 
Eq. (jnni), after the replacement (|H^ . to first order in as(/i) and find: 


Su;lf^\r,as) = l + ^Cp 

71 


4\n{N)Si{N)-2\n\N)-4+-TT^ + - - 2Si{N) Inr 


+ 0(« 


(33) 

This term correctly gives the large-iV pieces of the NLO cross section in Eq. (HID, but it goes 
beyond that by also reproducing all contributions \n(N)/N in Eq. ([HD, the latter arising from 
the expansion S'i(iV) = In(iV) + l/(2iV) + C>(l/iV^). 


The unpolarized resummed partonic cross section in Mellin-moment space is practically iden¬ 
tical to the transversely polarized one in Eq. dH, since the coefficients A{as) and Cq{r,as) are 
spin-independent and thus the same for the unpolarized and transversely polarized cases. A very 
small difference arises in the replacement in Eq. (El, in which for the unpolarized case one of 
course has to use the unpolarized LO splitting function. In addition, one should also take into 
account the singlet mixing in the evolution EH] which, however, is very small in the kinematic 
region we are interested in. 


3.2 Far-infrared resummed cross section 

We will see in the phenomenology section below that perturbative resummation as formulated 
so far in Eqs. dH or (j23|l predicts very large enhancements of the lowest order cross section, 
sometimes by orders of magnitude. There is good reason to believe that this enhancement is 
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only partly physical. The large corrections arise from a region where the integral in the exponent 
becomes sensitive to the behavior of the integrand at small values of kx- As long as Kqcd ^ 
M/N -C M, the use of NLL perturbation theory may be justihed, but when |A^| becomes so 
large that kr goes down to nonperturbative scales, we may well question the self-consistency of 
perturbation theory. We now turn to the question of how this “far-infrared” limit should be treated 
in resummed perturbation theory. Our discussion here will be brief, and we will only present one 
model that addresses the far-infrared limit. A more detailed study of this very interesting regime 
will be presented in a future publication. 


We seek a modihcation of the perturbative expression Eq. (P|) that excludes the region in which 
the absolute value of is less than some scale /xq > rriT,. We think of po as the scale beyond 
which the true mass spectrum of QCD replaces perturbation theory, regulating all soft and collinear 
singularities, so that should be thought of only as a lower limit for /tq. To implement this idea, 
we will adopt a modihed resummed hard scattering, which reproduces NLL logarithmic behavior 
in the moment variable N so long as M/N > fiQ, but “freezes” once M/N < / iq . If nothing else, 
this will test the importance of the region kr < ^qcd for the resummed cross section. If N were 
real and positive, we could simply replace the hrst exponent in (123 by 


rM 


dki 




■Aq{as{kT)) In 


Nkj 


(34) 


where 


p(a, h) = max(a, h) , 


(35) 


and where po then serves to cut off the lower logarithmic behavior. To provide an expression that 
can be continued to complex N^ we choose 


p(o,5) = (o'+ 


(36) 


with integer p. This simple form is consistent with the minimal expansion given above, and it 
also allows for a straightforward analysis of the ensuing branch cuts in the complex-A^ plane. 
Vanishing po corresponds to the standard minimal form (I23|) . For dehniteness, and for simplicity, 
we choose p = 2 in this paper. We will continue to use the expansions in Eq. (Oil) . but redehning 
A in Eq. by 

A = In (^1 + • (37) 

This form has the advantage that it generates only even power corrections in the ratio 
when this quantity is not too large [20] • We will investigate below the moderation of the pertur¬ 
bative increase provided by the cut-off /iq. 


3.3 Matching to the NLO cross section, and inverse Mellin transform 

As we have discussed above, the resummation is achieved in Mellin moment space. In order to 
obtain a resummed cross section in ^ space, one needs an inverse Mellin transform. This requires 
a prescription for dealing with the singularity in the perturbative strong coupling constant in 
Eqs. (IT21). ([22|) or in the NLL expansion, Eqs. (123,(EZj)- We will use the Minimal Prescription 
developed in Ref. EH], which relies on use of the NLL expanded form Eqs. (123,(1221), and on 
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choosing an appropriate contour in the complex-iV plane. In the standard minimal prescription, 
based on Eqs. (El and (El, this contour is chosen to lie to the left of the “Landau” singularities 
at A = 1/2 in the Mellin integrand, which are far to the right in the N plane. With the modihed 
variable A of Eq. El, however, branch cuts from A = 1/2 reside on the imaginary axis in the N- 
plane (and so do branch cuts that arise when the argument of the logarithm in Eq. (jHTjl vanishes). 
We again choose our inverse contour as 

drdcj) 

where Cmp is any positive number. All singularities are then to the left of the contour. We 
keep the Mellin contour parallel to the imaginary-A^ axis. The result dehned by the minimal 
prescription has the property that its perturbative expansion is an asymptotic series that has 
no factorial divergence and therefore no “built-in” power-like ambiguities. Power corrections may 
then be added, as phenomenologically required, and in a sense our cut-off prescription does exactly 
this. 


.CMP+ioo _^^^^(res),7V 




2711 


T 


d(f) 


(38) 


When performing the resummation, one of course wants to make full use of the available 
hxed-order cross section, which in our case is NLO {0{as))- Therefore, a matching to this cross 
section is appropriate, which may be achieved by expanding the resummed cross section to (P(as), 
subtracting the expanded result from the resummed one, and adding the full NLO cross section: 


drdcj) 




’‘CMP+ioo JT\T 






2711 


X 




0{as) 


+ 


drdcj) 


(39) 


In this way, NLO is taken into account in full, and the soft-gluon contributions beyond NLO are 
resummed to NLL. Any double-counting of perturbative orders is avoided. As we will see below, 
however, matching is almost academic in the present calculation since the NLO-expansion of the 
resummed cross section agrees to within 0.1% with the full NLO one for the kinematics relevant 
here. We also note that whenever we will use the form EH) with a non-vanishing cut-off po, we 

in Ei evaluated at po = 0. 


will perform the matching using the expansion as{fi)) 


0{a,) 


4 Phenomenological Results 


Starting from Eq. Ei, we are now ready to present some hrst resummed results at the hadronic 
level. This is not meant to be an exhaustive study; rather we should like to investigate the overall 
size and relevance of the resummation effects. For this reason, we only consider the cross section 
d'^a/dMdcj), integrated over all rapidities. This should be sufficient to study the main effects. In 
experiment one will eventually study rapidity distributions in order to better pin down the x- 
dependences of the transversity densities. Our resummation could be extended to this case using 
techniques developed in EH]- We also note that the charmonium resonances will dominate over a 
part of the spectrum that we will consider. For our case study we will just ignore this, but remind 
the reader that a complete treatment of the dilepton spectrum will eventually also require the 
incorporation of charmonium production and its resummation effects. 
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We will investigate pp collisions at four different energies. Each of these may eventually be 
realized at the GSI. The hrst two are in the hxed-target mode, relevant for the initial stage of pp 
physics at the GSI, when antiprotons will be scattered off proton targets at energies in the range 
15 GeV < Ep <25 GeV. For dehniteness, we will consider S' = 30 GeV^ and S' = 45 GeV^. This is 
also the regime in which the theoretical description is most challenging. The second regime is for 
an asymmetric collider, with Ep = 3.5 GeV and Ep = 15 GeV, corresponding to y/S = 14.5 GeV. 
Finally, we will consider a symmetric collider with Ep = Ep = 15 GeV {\/S = 30 GeV). In all 
calculations below, we choose the renormalization and factorization scales as p = M. 


4.1 Unpolarized cross section near partonic threshold 


We start by considering the unpolarized cross section. For some of the kinematic regions described 
above, Drell-Yan experiments at the GSI would enter uncharted territory, and it will be crucial 
to develop conhdence that the theoretical framework is understood. An important test would be 
a comparison to precise measurements of the unpolarized cross section. 

For our unpolarized calculations we will use the NLO (MS scheme) GRV parton distributions 
throughout jlU]- In the unpolarized case, we are in the fortunate situation that even the full 
NNLO corrections to the partonic Drell-Yan cross section are available BH, which we will in¬ 
corporate in our studies. These should in principle be used in conjunction with a set of NNLO 
parton distributions, which became available recently 1321 after the computation of the three-loop 
evolution kernels jl21- The main purpose of our present studies, however, is to see how well the 
soft-gluon terms we are resumming reproduce also the NNLO corrections to the cross section. For 
this it is sufficient to stick to our use of NLO parton densities. 


Figure n shows the effects of the higher orders generated by resummation, for the hxed-target 
cases with S' = 30 GeV^ and S = 45 GeV^. We dehne a resummed “iL-factor” as the ratio of the 
resummed (matched) cross section to the LO cross section. 




(40) 


which is shown by the solid line in Fig. ^ As can be seen, is very large, meaning that resum¬ 

mation results in a dramatic enhancement over LO, sometimes by over two orders of magnitude. 
It is then interesting to see how this enhancement builds up order by order in the resummed cross 
section. We expand the matched resummed formula to NLO and beyond and dehne the (here, 
not matched) “soft-gluon iL-factors” 


da^^^'>/dMd(j) ’ 


(41) 


which for n = 1 , 2 ,... give the additional enhancement due to the 0{a'^) terms in the resummed 
formula. Formally, = 1, while up to the ehects of matching at NLO. The results 

for iLi’2>3,4,6,8 in Fig. ^ One can see that there are very large contributions even 

beyond NNLO, in particular at the higher M. Glearly, the full resummation given by the solid 
line receives contributions from high orders. The symbols in Fig. C] show the associated K factors 
for the exact NLO and NNLO calculations, respectively. One can see that these agree extremely 
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Figure 1: “iF-factors” relative to LO as defined in Eqs. (gni) and dUD for the Drell-Yan cross 
section in fixed-target pp collisions at S' = 30 GeV^ (left) and S = 45 GeV^ (right), as functions 
of lepton pair invariant mass M. The symbols denote the results for the exact NLO and NNLO 
calculations. 




Figure 2: Same as Fig.[Tl but for pp collider options with a/S = 14.5 GeV (left) and a/S = 30 GeV 
(right). 


well with the 0{as) and O(a^) expansions of the resummed cross section. In fact, the agreement 
between the full NLO result and the 0{as) expansion of the resummed cross section is better than 
0.1% over the whole range in M shown. Thus the matched resummed iF-factor of Eq. (1401) 

is also numerically very close to iF°° as defined in dUD- We note that the replacement in Eq. (El 
helps somewhat in this comparison, in particular at NNLO it leads to a relative improvement of 
a few per cent. From the comparison we may conclude that the terms that we resum to all orders 
strongly dominate the cross section. 

In Fig. H we show similar results for the two collider modes at VS = 14.5 and 30 GeV. One 
can see that for a hxed Drell-Yan mass M the corrections become much smaller, since one is much 
further away from partonic threshold than for the hxed-target cases. Also, the convergence of the 
perturbative series occurs somewhat more rapidly. As before, the agreement between the 0(as) 
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Figure 3: Unpolarized cross sections dPa/dM at S' = 30 GeV^ (left) and S' = 45 GeV^ (right) at 
LO, NLO, NNLO, and NLL resummed, as functions of lepton pair invariant mass M. 



Figure 4: Same as Fig.lHl but for pp collider options with a/S = 14.5 GeV (left) and a/S = 30 GeV 
(right). 


and 0{al) expansions of the resummed cross section and the exact NLO and NNLO calculations 
is very good, demonstrating the relevance of the resummed result. 

In Figs. Elandm we show the actual unpolarized cross sections M^da/dM corresponding to the 
results shown in Figs. [Hand 121 at (exact) hxed orders (LO, NLO, NNLO) in perturbation theory, 
and for the NLL resummed case. Here we have integrated over all azimuthal angles 0'. 

One may wonder whether any sign of the need for large corrections beyond NLO can be found 
in previous Drell-Yan data jH]. The measurements at the lowest '/S we are aware of are from 
the GERN WA39 |12] experiment and were made in vr^-Tungsten scattering. The pion energy 
was Ett = 39.5 GeV, higher than what is considered for the hxed-target mode at the GSI, but 
still quite far below the collider energies. Fig. El shows the data, along with our results at LO, 
NLO, NNLO and NLL-resummed. We are using the parton distributions for the pion of jlH]- It is 

^Note that the normalization factor in Eq. © becomes 20^6^193 in the unpolarized case. 
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Figure 5: Comparison of LO, NLO, NNLO and resummed results to the WA39 data jlH] for 
Drell-Yan dimuon production in scattering. 


hard to draw definite conclusions from the comparison in Fig. 0 partly because the experimental 
uncertainties are rather large, and also because the pion parton densities are not known accurately. 
Nevertheless, it is interesting to observe that the data points at the highest dimuon mass M are 
quite consistent with large perturbative resummation effects. 

Given the large, not to say huge, size of some of the enhancements, we now turn to the 
same ratios computed using the exponent (I34|l with the modified lower limit on the integral, 
which regulates its infrared behavior. For simplicity, we will only show results for relatively small 
values of the cut-off scale /io in Of course, different choices give different results, but we 

should think of /xq ns a kind of factorization scale, separating perturbative contributions from 
nonperturbative. Thus changes in /xq would be compensated at least in part by changes in a 
nonperturbative function. On the other hand, a very strong sensitivity to /xq can reasonably be 
interpreted as indicating that perturbative resummation alone cannot give a reliable estimate for 
the cross section. Our interest here, therefore, is primarily to illustrate the modihcation of the 
perturbative sector, which we do by choosing /xq = 0.3 GeV and 0.4 GeV. Results for the “iF- 
factor” with these values of /xq are shown in Figs. El and [71 compared to the same NLO, NNLO 
and resummed factors shown above at the relevant energies and pair masses M. We hope to 
study the /xq dependence more extensively elsewhere, in connection with possible nonperturbative 
corrections. 

The ratios of the new resummed but infrared-regulated cross sections to the LO one show a 
smoother increase than the pure minimal resummed cross sections. This difference is particularly 
marked at the lower center-of-mass energies in Fig. El with only a modest enhancement over 
NNLO remaining at /xq = 0.3 GeV, and even lower at 0.4 GeV. We note that the NLO and NNLO 
expansions of the resummed cross section turn out to be much less affected by the cut-off /xq than 
the full resummed cross section. At the higher energies of Fig. [71 the regulated resummed curves 
follow the unregulated curves far above NNLO, with much reduced sensitivity to /xq. We interpret 
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Figure 6: “iF-factors” relative to LO as in Fig. [H at S = 30 GeV^ (left) and 5 = 45 GeV^ 
(right). The dashed lines show the effect of a lower cut-off /iq = 300 MeV for the /cT-integral in 
the exponent, and the dot-dashed a cut-off of /xq = 400 MeV for comparison. The symbols denote 
the results for the exact NLO and NNLO calculations. 




Figure 7: Same as Fig.lHl but for pp collider options with y/S = 14.5 GeV (left) and y/S = 30 GeV 
(right). 


these results to indicate a strong sensitivity to nonperturbative dynamics at the lower energies, 
and much less at the higher. We therefore take the predictions of large enhancements due to 
high order perturbation theory more seriously in the latter case, even while keeping in mind the 
WA39 measurements of Fig. El above, which suggest that large corrections, perturbative or not, 
should not be ruled out a priori. Data over the entire kinematic regime considered here would 
certainly shed a unique light on the transition between long- and short-distance effects in hadronic 
scattering. 
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4.2 Spin asymmetry Aj'j' 


Before we can perform numerical studies of A^t we need to make a model for the transversity 
densities in the valence region. Here, guidance is provided by the Softer inequality 


2 < q(x,Q^) + Aq(x,Q^) , 


(42) 


which gives an upper bound for each 6q. Following daini we utilize this inequality by saturating 
the bound at some low input scale Qo ~ 0.6 GeV using the NLO GRV [10] and GRSV (“standard 


scenario”) jlH] densities q{x, Ql) and Aq{x, Ql), respectively. For Q > Qo the transversity densities 
Sq{x, Q^) are then obtained by solving the evolution equations with the NLO pTfl evolution 
kernels. We refer the reader to OEI for more details on our model distributions. 

We will now investigate to what extent the large perturbative corrections we found for the un¬ 
polarized Drell-Yan cross section cancel in the spin asymmetry A^t, Eq- ©■ We have mentioned 
earlier that the resummation factors (HI for the qq cross section are the same in the unpolarized 
and transversely polarized cases if both are treated in the same factorization scheme. Resum¬ 
mation effects would therefore cancel if the spin asymmetry were in Mellin-moment space. The 
convolution with the parton distributions and the inverse Mellin transform will affect the cancel¬ 
lation somewhat, but one still expects the spin asymmetry to be very robust. Indeed, the results 
for Att at LO, NLO, and resummed to NLL (with and without the cut-off /iq), shown in Figs. |H1 
and El for the four energies that we consider, confirm this. We show 


Aj'j' — 


d5a/dMdcj) 


(43) 


da/dMdcj) 

as a function of M. For simplicity we set 0 = 0 here; extension to other cj is straightforward by 
taking into account the cos(20)-dependence displayed in Eq. (jHl). 

Aside from a slight deficit in the curves at lower M (and slight excess at higher M) with 
unmodihed minimal resummation at the lower S, all curves, including the LO, NLO and regulated 
resummed asymmetries all he within a few percent of each other. We note that to NLO this 
robustness of Att was also found for pp collisions at higher energies lag. 

We can shed some light on why the asymmetry is modestly but significantly shifted for the 
unregulated perturbative resummed cross section with respect to the other cases shown. First, we 
observe that smaller momentum fractions Xa^t, in the factorized cross section are associated with 
the cross section at smaller pair mass M, where the spin asymmetry is slightly smaller. Therefore, 
any effect that tends to drive momentum fractions particularly close to their minimum values, at 
z = 1, will tend to decrease the asymmetry. We note, however, that because the valence quark 
distributions are decreasing functions of the x’s, and because they are the same for both hadrons, 
we might anticipate that the average values (xa) = (xf,) are not far from their symmetric values, 
Xa = Xb = a/t at partonic threshold. In fact, this turns out to be a surprisingly accurate estimate, 
even at lowest order, as can be readily verihed from Fig. E3 where we plot (x) for LO, NLO, and 
the unregulated (po = 0) and regulated (po = 300 MeV) resummed cross sections. In this hgure 
(x) is found by performing the integral for the unpolarized cross section with an extra factor Xa 
in the integrand, and by dividing by the cross section itself. 

Beyond lowest order, the hard-scattering cross section is further enhanced at partonic threshold 
{z = 1), and we would expect that large perturbative enhancements, such as those associated 
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Figure 8: Spin asymmetry Att{(P = 0) LO, NLO and for the NLL resummed case, at S' = 
30 GeV^ (left) and S = 45 GeV^ (right). The dash-dotted lines show the effect of a lower cut-off 
fio = 300 MeV for the fc-r-integral in the exponent. 




Figure 9: Same as Fig.|Hl but for pp collider options with y/S = 14.5 GeV (left) and y/S = 30 GeV 
(right). 


with the plus distribution of Eq. would force the average partonic center-of-mass energy still 
closer to threshold, and hence reduce the asymmetry even further. We should observe, however, 
that the hard scattering function is not a positive-dehnite cross section, but rather a sum of 
plus distributions, given at hrst nontrivial order in Eq. ©• At NLO, for example, the positive 
contribution at z = 1 is from a delta function associated with virtual corrections, while the real- 
gluon contribution, ln(l — z)/{I — z), is actually negative, due to the subtraction of collinear 
divergences in the calculation of the hard scattering |221I21]- We therefore cannot interpret {xa,b) 
as averages in the usual sense. In any case, we do see a more signihcant decrease in (x), computed 
above, for the unregulated resummed cross section than for fixed order. In fact, the values derived 
in this manner are below which would be the lower limit for a positive-dehnite hard scattering 
function, and at the highest M, even below r, which is the lower limit of the integration range 
for the x’s. The caveat against a literal interpretation of {x) notwithstanding, it is reasonable to 
interpret the modest decrease in the asymmetry for the unregulated resummed cross section as 
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Figure 10: Average values {xa) = {xh) = (x) for the parton momentum fractions at S' = 30 GeV^ as 
functions of M, at LO, NLO, for purely perturbative resummed, and for regulated resummed with 
/To = 300 MeV. For comparison, the lowest curve shows the corresponding values of r = M^/S. 
A curve for a/t - a straight line in this plot - would be practically indistinguishable from the LO 
one, although remaining slightly below it. 


resulting from a hard-scattering function that is exceptionally peaked near z = 1 in this case. 


4.3 Use of DIS parton distribution functions 


We recall that throughout this study, we have used the NLO parton distribution functions of 
Ref. |in]. Especially given the large effects we have found with resummation, we should revisit 
the use of “un-resummed” parton distributions in our study. The momentum fractions probed in 
the parton distributions become very large in the hxed-target regime, as shown by Fig. One 
may wonder here to what extent the parton distributions themselves should include resummation 
effects. The densities we use have been determined mostly from an analysis of data from deeply- 
inelastic scattering (DIS), in which however no resummation of large-W logarithms was included. 
It is known that soft-gluon resummation effects in DIS are rather unimportant, except at moderate 
photon virtuality and very large x ISIIE21 Nonetheless, we will give a rough estimate of the 
quantitative effect on the Drell-Yan cross section that might occur if one used parton densities 
determined from an analysis of the DIS data including resummation. We follow ISH to determine a 
model set of “MS-resummed” valence distributions gNres Mellin-moment space) by demanding 
that their contributions to the DIS structure function F 2 match those of the corresponding NLO 
densities at a hxed scale Q. This is ensured by “rescaling” the parton densities: 


^N.res 


(Q-) = 


Gr(iV,Q2) 


(44) 
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where and Cf" are the perturbative NLO and NLL resummed quark coefficient functions for 
F 2 , respectively, which may be found in for example. We choose Q in Eq. (I44j) fairly large, so 
that it is in a region where resummation effects are expected to be small and NLO to yield a good 
description of F 2 . For illustration, we use two different values, = 25 GeV^ and = 100 GeV^. 
The ansatz (m then represents an estimate of the likely change in the parton distributions from 
resummation in DIS, and the parton densities gNres used for calculating the Drell-Yan 

cross section. The result is shown in Fig. where we repeat the resummed cross section from 
Fig.El and display the effects on the cross section for the two choices of Q. A moderate decrease 
of the resummed cross section is found. We note that this effect becomes smaller at the collider 
energies, for a given M. 



= 30 GeV^ due to “rescaling” 
result for the resummed cross 
and 10 GeV, respectively, in 


5 Conclusions 


Figure 11: Effect on the unpolarized Drell-Yan cross section at S 
of the parton distributions as in Eq. dm). The solid line shows the 
section as in Fig. El the dotted and dashed lines are for Q = 5 

Eq. (jUj). 


We have verihed that perturbative corrections associated with partonic threshold are large for 
dilepton production in the kinematic region being considered for proton-antiproton collisions at 
GSI. The close agreement between the fixed-order expansions of threshold-resummed cross sec¬ 
tions and exact fixed order calculations suggests that the resummation is physically relevant, and 
signihcantly enhances the unpolarized Drell-Yan cross sections beyond hxed orders. At the same 
time, the transverse-spin asymmetries whose measurement is suggested for these experiments, are 
remarkably insensitive to shifts in the overall normalization. In summary, perturbative corrections 
appear to make the cross sections larger independently of spin. They would therefore make easier 
the study of spin asymmetries, and ultimately transversity distributions. 
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We have also shown that at the lower energies considered, the resummed cross section decreases 
markedly when an infrared cut-off is used to regulate contributions from soft gluon emission in 
the far infrared region. We have shown how to incorporate such a cut-off as a generalization of 
minimal resummation. At the lower energies, a strong sensitivity to the cut-off suggests that the 
large enhancements found with unregulated perturbative resummation arise from an unwarranted 
extension of perturbation theory into the soft region At the same time, it is interesting to note that 
quite substantial A'-factors survive almost unchanged by infrared regulation at the higher energies 
considered. This suggests that the measurement of the unpolarized dilepton cross section will shed 
light on the relationship between hxed orders, perturbative resummation and nonperturbative 
dynamics in hadronic scattering. 
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